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condition  EP3).  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent  CiC/MS 
data. 

Fig.  9b  Minor  species  (<  1 0  ppm)  for  pyrolysis  of  ethyl  methyl  sulfide  at  740  X  (operating 
condition  EP3).  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent  GC/MS 
data. 

Fig.  10a  Major  species  <>!0  ppm)  for  oxidation  of  diethyl  sulfide  at  630  X  with  low  02  loading 
(operating  condition  DOl).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 
GC/MS  data. 

Fig.  10b  Minor  species  (<10  ppm)  for  oxidation  of  diethyl  sulfide  at  630  °C  with  low  0:  loading 
(operating  condition  DOl).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  FTIR  data  only. 

Fig.  1  la  Major  species  (>10  ppm)  for  oxidation  of  diethyl  sulfide  at  670  X’  with  low  02  loading 
(operating  condition  D<)2).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 
GC/MS  data. 

Fig.  I  lb  Minor  species  (<10  ppm)  for  oxidation  of  diethyl  sulfide  at  670  X  with  low  02  loading 
(operating  condition  D<)2).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  FTIR  data  only. 

Fig.  12a  Major  species  <>!0  ppm)  for  oxidation  of  diethyl  sulfide  at  700  "C  with  low  02  loading 
(operating  condition  D<)3).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 
GC/MS  data. 

Fig.  12b  Minor  species  (<10  ppm)  for  oxidation  of  diethyl  sulfide  at  700  °C  with  low  ()2  loading 
(operating  condition  D()3).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  FTIR  data  only. 

Fig.  13a  Major  species  <>!0  ppm)  for  oxidation  of  diethyl  sulfide  at  740  "C  with  low  <)2  loading 
(operating  condition  D<)4).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 
GC/MS  data. 

Fig.  14a  Major  species  (>10  ppm)  for  oxidation  of  diethyl  sulfide  at  700  X  with  high  O.  loading 
(operating  condition  DOS).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 
GC/MS  data. 

Fig.  14b  Minor  species  (<I0  ppm)  for  oxidation  of  diethyl  sulfide  at  700  "C  with  high  02  loading 
(operating  condition  DOS).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  FTIR  data  only. 

Fig.  15a  Major  species  (>  1 0  ppm)  for  oxidation  of  diethyl  sulfide  at  740  X  with  medium  02 
loading  (operating  condition  D<)6).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols 
represent  GC/MS  data. 

Fig.  15b  Minor  species  (<  1 0  ppm)  for  oxidation  of  diethyl  sulfide  at  740  X'  with  medium  02 
loading  (operating  condition  D()6).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  FTIR  data  only. 
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Fig.  16a  Major  species  (>  1 0  ppm)  for  oxidation  of  diethyl  sulfide  at  740  *C  with  high  Cfc  loading 
(operating  condition  D07).  Lines  represent  calculations  with  mechanism  version  2;  symbols 
represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 
GC7MS  data. 

Fig.  16b  Minor  species  (<1()  ppm  measured)  for  oxidation  of  diethyl  sulfide  at  740  °C  with  high 
O:  loading  (operating  condition  D()7).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  FTIR  data  only. 

Fig.  17a  Major  species  (>10  ppm)  for  oxidation  of  ethyl  methyl  sulfide  at  700  °C  with  high  Cfc 
loading  (operating  condition  EOl).  FTIR  data  only 

Fig.  17b  Minor  species  (<10  ppm)  for  oxidation  of  ethyl  methyl  sulfide  at  700  ”C  with  high  O- 
loading  (operating  condition  EOl).  FTIR  data  only 

Fig.  18a  Major  species  (>10  ppm)  for  oxidation  of  ethyl  methyl  sulfide  at  740  "C  (operating 
condition  E<)2).  FTIR  data  only. 

Fig.  18b  Minor  species  (<10  ppm)  for  oxidation  of  ethyl  methyl  sulfide  at  740  "C  (operating 
condition  E<)2).  FTIR  data  only. 

Fig.  19  Element  balances  for  diethyl  sulfide  pyrolysis  data  sets  (operating  conditions  DPI.  DP2. 
DP3,  and  DP4.) 

Fig.  20  Major  reaction  pathways  for  diethyl  sulfide  pyrolysis.  Dotted  lines  indicate  radical  attack 
reactions.  Thicker  lines  represent  reactions  with  larger  rates. 

Fig.  21  Major  reaction  pathways  for  diethyl  sulfide  oxidation.  Dotted  lines  indicate  radical  attack 
reactions.  Thicker  lines  represent  reactions  with  larger  rates. 
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III.  PROBLEM  STUDIED 


We  report  the  results  of  a  study  investigating  the  destruction  chemistry  of  organosulfur 
compounds  under  pyrolytic  and  oxidative  conditions.  Sulfur  compounds  arc  common  in 
combustion  systems  because  of  the  presence  of  sulfur  impurities  in  fuels  such  as  coal  and 
oil.  However,  the  understanding  of  sulfur  combustion  chemistry  is  largely  limited  to  the 
pathways  leading  to  sulfur  oxides.  We  address  the  early  stages  of  destruction  for  alkyl 
sulfides,  which  are  surrogates  for  the  chemical  warfare  agents  H.  HD.  and  HT.  We  report 
computational  chemistry  and  experimental  flow  reactor  results  contributing  to  the 
development  and  testing  of  a  detailed  chemical  kinetic  mechanism  for  the  combustion 
and  pyrolysis  of  diethyl  sulfide.  This  mechanism  can  serve  as  a  starting  point  for 
developing  a  mechanism  for  the  related  chemical  warfare  agents,  which  would  be  of 
value  for  predicting  incinerator  performance  including  off-design  conditions.  A  similar 
effort  with  organophosphorus  compounds  has  led  to  a  detailed  kinetic  model  of  chemical 
warfare  agent  combustion,  and  its  use  in  a  model  of  an  incinerator  [  1 J. 

IV.  IMPORTANT  RESULTS 


This  section  describes  our  major  results.  First  we  briefly  summarize  the  computational 
chemistry  methods  and  results  for  thermochemical  properties,  reaction  pathways,  and 
kinetic  parameters.  Then  we  discuss  the  two  versions  of  the  sulfide  oxidation/pyrolysis 
mechanism  developed  during  the  project.  Then  we  discuss  the  flow  reactor  experimental 
apparatus  and  method  as  well  as  the  method  for  performing  chemical  kinetic  simulations 
of  the  flow  reactor  experiments.  Finally,  experimental  results  and  mechanism  predictions 
arc  presented. 

IN'. A.  Computational  Chemistry  and  Mechanism  Development  Results 

IV.A.l.  Calculations  of  thennochemical  properties 


We  have  computed  thermochemical  properties  for  several  classes  of  sulfur-containing 
species.  The  following  method  was  used  for  the  calculations:  Isodesmic  work  reactions 
with  various  levels  of  theory  have  been  used  to  calculate  enthalpies  of  formation 
(AHr°29s)  of  stable  species  and  radicals.  Entropy  (S)  and  Heat  capacity  (CP(T))  arc 
determined  using  geometric  parameters,  vibrational  frequencies  and  moments  of  inertia. 
Internal  rotor  potentials  are  calculated  at  the  same  level  of  theory. 

Table  1  shows  the  classes  of  compounds  for  which  thermochemical  properties  have  been 
calculated,  along  with  the  theory  level  used.  The  thennochemical  data  resulting  from 
these  calculations  arc  reported  in  our  previous  progress  reports  as  well  as  in  our 
publications  and  presentations.  Thermochemical  properties  calculated  as  parts  of  the 
kinetic  investigations  described  in  the  next  section  arc  not  included  in  Table  1. 
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Table  I.  (.'lasses  of  compounds  for  w  hich  thermochcmical  properties  have  been  obtained 
through  computational  chemistry. 


Class  of  compounds 

Example  species 

Level  of  theory 

Sulfcnic  acids  and  esters 
and  associated  radicals 

CHj-S-OH,  CH3-S-OCH3, 
CH2-S-OH 

CBS-QB3,  G3MP 

Peroxy  and  hydroperoxy 
radicals  associated  with 
disulfide  molecules 

HSS-CH:-OOH,  CHs-SS- 
CH2-OOH,  -SS-CFL-OOH. 
HSS-CFL-OO 

B3  LY  P/6-3  l(d,p), 

CBS-QB3 

Methyl  sulfinic  acid  and 
ester  and  associated  radicals 

CH3-S(=0)-0H, 

CH3-S(=0)-0CH3, 

CH3-S(=0)-0* 

B3  L  Y  P/6-3 1  G(2d.p), 

B3LYP/6-31G(2d.dp), 

CBS-QB3 

tthencthiol  and  related 

radicals 

CH2=CH-S-H.  CH2=CH-S* 

CBS-QB3.  G3MP 

IV.A.2.  Calculations  of  reaction  pathways  and  kinetics 

For  a  number  of  reactions  of  interest  in  the  hydrogcn/carbon/sulfur/oxygen  system,  an 
initial  reaction  forms  a  chemically  activated  adduct  to  which  numerous  reaction  pathways 
arc  available.  This  set  of  reaction  pathways  has  been  investigated  computationally  for  a 
number  of  initial  reactions,  listed  in  Table  2. 


Table  2.  Multichannel  reactions  for  which  energy  levels  and/or  kinetic  parameters  have 
been  obtained  through  computational  chemistry. 


Initial  reaction 

Level  of  theory 

HS\+CO 

CH,S-  +CO 

CHjCHjS*  +CO 

CH; = CH-  CH:S*  +CO 

CBS-QB3  multilevel  composite  method  for  smaller 
radicals;  density  functional  theory  at  B3LYP/6- 
31G(d.p)  level. 

CHsSCH +  'O: 

CBS-QB3  and  B3LYP/6-31  IG(d,p)) 

CHiS*  +  '(); 

Multiple  methods  including  B3LYP/6- 
31 1++G(d.p),  B3LYP/6-31 1++G(3df,2p), 
CCSD(T)/6-3 1  lG(d.p)  //MP2/6-3  lG(d,p), 

B3  P86/6-3 1 1  G(2d,2p)//B3  P86/6-3 1  G(d),  and 

B3PW9 1/6-3 1 1++G  (3df.2p).  Tlie  CBS-QB3  and 
G3MP2  composite  methods  are  recommended.  G2 
and  G3  arc  also  applied  when  necessary. 

CBS-QB3  and  density  functional  theory 

S=CmH+  !02 

( oxidation  of  thioformyl 
radical ) 

CBS-QB3,  G3MP2,  and  G2  for  enthalpies  of 
reaction;  B3LYP/6-3 1  lG(d.p)  for  S  and  Cp 
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The  procedure  for  obtaining  reaction  pathways  and  kinetics  is  the  following:  First 
thermochcmical  properties  are  calculated  for  the  adducts  and  the  products  of  the  reaction 
paths,  as  well  as  for  the  transition  states  for  all  reactions.  The  method  for  these 
calculations  was  similar  to  that  described  in  the  previous  section.  These  results  can  be 
displayed  in  a  diagram  showing  computed  potential  energy  levels  for  each  set  of  products 
or  transition  state.  This  potential  energy  level  diagram  can  be  interpreted  qualitatively, 
to  sec  which  channels  have  relatively  low  barriers,  and  thus  arc  likely  to  be  most 
important.  A  representative  energy  level  diagram  is  shown  below,  as  Figure  1. 


kcnl  t  mol  it-ii 


Figure  I:  Potential  Energy  Diagram  at  298  K  (in  kcaFmol)  for  the  oxidation  of  the 
thioformyl  radical.  All  harriers  are  calculated  at  C3.MP2  level,  except  those  labeled  with  * 
and  **.  which  are  calculated  at  CBS-QB3  and  the  average  of  several  levels,  respectively. 

In  the  cases  where  more  quantitative  kinetic  data  was  obtained,  it  was  done  so  as  follows: 
The  rate  constants  for  the  elementary  reactions  were  calculated  from  the  computed 
frequencies,  structures  and  energies  using  canonical  transition  state  theory  with  rate 
constant  values  (T)  Fit  to  the  three-parameter  modified  Arrhenius  equation,  k  =  A  x  T°  x 
cxp(-Ea'RT).  Temperature  and  pressure  dependant  rate  constants  were  calculated  with 
multi  channel,  multi  frequency  QRRK  theory'  for  k{E)  and  Master  equation  analysis  for 
fall-off  [2].  A  set  of  rate  constants  calculated  from  the  energy  levels  shown  in  Figure  1  is 
shown  below  in  Figure  2.  Results  of  this  analysis  and  others  are  presented  in  our  reports 
and  publications. 
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Figure  2:  Calculated  rate  constants  ( l«y  k)  vs  temperature  (K)  at  I  atm  for  reactions 
inxolved  in  the  oxidation  of  the  thiformyl  radical. 

71  '..-1.3.  Development  of  a  combustion  chemistn'  mechanism  for  diethyl  sulfide 

A  pyrolysis  and  oxidation  mechanism  has  been  developed  for  diethyl  sulfide,  along  with 
species  thcrmochcmical  parameters  (enthalpy,  entropy,  and  heat  capacity)  for  molecule, 
radical  intermediates,  and  transition  state  structures.  The  pyrolysis  submechanism  of 
Version  I  of  the  mechanism  consists  of  1 70  elementary  reactions  and  70  species  that  are 
components  of  approximately  25  chemical  activation  or  multi-channel  dissociation 
reaction  systems.  This  submechanism  is  described  in  detail  in  Zheng  et  al.[3].  Version  2 
is  a  similar  mechanism  with  computational  chemistry’  used  at  high  levels  to  calculate  rate 
constants  on  important  reactions  of  the  reacting  sulfides  and  thiols.  Version  2  includes 
rate  constants  determined  at  the  CBS-QB3  and  G3MP2  levels  of  theory  for  abstraction, 
elimination  (beta  scission  from  radical  species  and  molecular  elimination  from  both 
radicals  and  singlet  molecules)  and  addition  reactions  related  to  the  parent  diethyl  and 
dimethyl  sulfides  and  ethyl  and  methyl  thiols.  Reactions  for  abstraction  and  elimination 
of  hydrogen  atoms  from  the  sulfide  and  thiol  moieties  include  tunneling  corrections  using 
the  Wigncr  method,  plus  degeneracy  and  symmetry  contributions. 
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IV. B.  Experimental  Results  and  Comparison  of  Mechanisms  lo  Experiments 


IV.B.I.  Experimental  method 

We  have  used  a  turbulent  flow  reactor  with  extractive  sampling  to  examine  the 
destruction  of  two  mustard  simulants  under  both  pyrolytic  and  oxidative  conditions.  The 
flow  reactor  is  shown  in  Fig.  3  and  is  described  briefly  here;  more  details  are  given  in  a 
recent  publication  [3J.  The  flow  reactor  consisted  of  a  main  flow  and  a  secondary  flow. 
Nitrogen  was  metered  and  passed  through  a  retort  inside  a  furnace.  The  main  flow  went 
through  a  heated  flow  development  section,  followed  by  an  injection  section  in  which  the 
secondary  flow  (simulant  +  N’  and  O;  or  CO)  was  added  via  four  turbulent  jets.  Next 
came  the  reaction  section,  a  quartz  tube  with  45  mm  ID.  48  mm  OD,  and  100  cm  length, 
held  inside  a  larger  heated  stainless  steel  pipe.  Four  ports  on  the  stainless  steel  pipe  were 
aligned  with  the  holes  in  the  quartz  tube,  allowing  the  insertion  of  probes  for  centerline 
gas  sampling.  Temperatures  measured  along  the  flow  reactor  axis  are  shown  in  Fig.  4. 


=1 


Figure  3.  Schematic  of  Flow  Reactor 
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Figure  4.  Reaction  section  axial  temperature  profiles.  The  four  different  symbols  represent 
different  operating  conditions  with  nominal  temperatures  of  630.  670.  700.  and  740  C. 

Lower  temperatures  occur  at  short  distances  from  the  injectors  as  the  cooler  secondary  flow 
mixes  into  the  main  flow. 

Samples  were  withdrawn  from  the  reactor  via  a  quartz  probe.  In  order  to  quench 
reactions,  the  sample  flow  was  diluted  to  approximately  50%  of  its  original  concentration 
with  a  flow  of  nitrogen  inside  the  quartz  probe.  A  simple  adiabatic  mixing  calculation 
indicates  that  the  dilution  should  cool  the  sample  flow  to  temperatures  at  which  no 
composition  changes  in  the  product  stream  are  predicted  to  occur.  Downstream  of  the 
probe,  the  diluted  sample  was  drawn  through  Teflon®  (Dupont)  tubing  cither  into  the 
long-path  gas  cell  of  a  FTIR  (Nicolct  6700)  or  into  the  gas  sampling  loop  of  a  GC/MS 
(Thermo  Focus  GCDSQll  MS)  for  sample  analysis.  For  details  of  the  sampling  and 
analysis  methods,  sec  Zheng  ct  al.  [3].  The  FTIR  was  used  mainly  for  light 
hydrocarbons,  while  the  GC/MS  yielded  information  about  the  parent  compound  as  well 
as  minor  sulfur-containing  products.  Typical  sensitivities  of  the  two  instruments  range 
from  50  ppb  to  1.5  ppm. 

The  sample  species  concentrations  depend  on  the  sample  dilution  and  on  the  rate  of 
mixing  of  main  and  secondary  flows,  as  well  as  on  the  rates  of  chemical  reactions 
occurring  in  the  flow  reactor.  For  the  pyrolysis  experiments,  CO  was  shown  to  be  inert, 
and  was  used  as  a  tracer  to  adjust  the  measured  species  concentrations  to  account  for 
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mixing  and  dilution.  The  measured  mole  fraelion  CO  versus  distance  is  shown  in  Fig.  3 
for  the  different  temperature  conditions.  For  the  oxidation  experiments,  CO  was  not  a 
viable  tracer,  as  it  could  react  with  oxygen.  The  mixing/dilution  information  obtained 
from  CO  tracer  during  the  pyrolysis  experiments  was  consistent  from  run  to  run,  and  was 
applied  to  the  oxidation  data  sets. 

IV.B.2.  Computational  method 

The  How  reactor  was  simulated  as  a  plug  How  reactor  in  CHEMKIN  [4],  with  centerline 
measured  temperatures  as  inputs  to  the  calculation,  and  with  no  wall  reactions. 
Experiments  were  performed  with  either  diethyl  sulfide  or  ethyl  methyl  sulfide  as  the 
parent  compound.  (See  operating  conditions  below.)  Only  the  diethyl  sulfide 
experiments  were  simulated,  as  the  kinetic  mechanism  did  not  include  the  initial  stages  of 
the  decomposition  of  ethyl  methyl  sulfide.  Two  different  kinetic  mechanisms  are  used, 
as  described  in  Section  IV. A. 3  above.  Mechanism  Version  1 .  described  in  a  previous 
publication  [2J,  was  used  to  simulate  the  diethyl  sulfide  pyrolysis  experiments. 
Mechanism  Version  2,  which  differs  from  Mechanism  Version  1  only  through  the  use  of 
a  higher  level  of  theory  in  the  calculation  of  some  rate  constants  and  thermochemical 
parameters,  was  used  to  simulate  the  diethyl  sulfide  oxidation  experiments. 
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II  B.3.  Operating  conditions 

Tables  3  through  6  show  the  flow  reactor  operating  conditions  for  the  four  experimental 
data  sets  obtained  (total  of  12  operating  conditions  with  two  different  chemical  warfare 
agent  simulants).  These  data  sets  arc:  pyrolysis  of  diethyl  sulfide  (Table  3)  and  of  ethyl 
methyl  sulfide  (Table  4),  and  oxidation  of  diethyl  sulfide  (Table  5)  and  of  ethyl  methyl 
sulfide  (Table  6).  The  temperature  profiles  are  as  shown  in  Figure  4.  In  tables  5  and  6, 
the  equivalence  ratio  is  based  on  the  conversion  of  sulfur  to  SO2,  carbon  to  CO2,  and 
hydrogen  to  H^O  for  complete  combustion 


Table  3.  Flow  reactor  operating  conditions  for  diethyl  sulfide  pyrolysis  experiments.  In 
these  experiments  the  secondary  flow  consisted  of  N2,  CO,  and  diethyl  sulfide. 


Condition 

Nominal 

Temperature 

rc> 

Flowrate 
of  Main  Nj 
flow 
(slpm) 

Flowrate  of 
Secondary 

N-  flow 
(slpm) 

Well- 

mixed 

CO 

loading 

in 

reaction 

section 

(ppm) 

Well-mixed 
initial  diethyl 
sulfide  loading 
(ppm) 

Reynolds 

number 

DPI 

630 

330 

14.1 

68.6 

150 

5022 

DP2 

670 

341 

15.3 

66.3 

150 

5022 

DP3 

700 

347 

16.0 

65.1 

150 

5022 

DP4 

740 

350 

16.8 

64.4 

150 

5022 

Table  4.  Flow  reactor  operating  conditions  for  ethyl  methyl  sulfide  pyrolysis  experiments. 
In  these  experiments  the  secondary  flow  consisted  of  N;,  CO.  and  ethyl  methyl  sulfide. 


Condition 

Nominal 

Temperature 

("C> 

Flowrate 
of  Main 

N;  HOW 
(slpm) 

Flowrate  of 
Secondary 

N.  flow 
(slpm) 

Well- 

mixed 

CO 

loading 

in 

reaction 

section 

(ppm) 

Well-mixed 
initial  ethyl 
methyl  sulfide 
loading  (ppm) 

Reynolds 

number 

EP1 

630 

330 

14.1 

68.6 

150 

5020 

EP2 

700 

347 

16.0 

65.1 

150 

5060 

EP3 

740 

350 

16.8 

64.4 

150 

49% 
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Tabic  5.  Flow  reactor  operating  conditions  for  diethyl  sulfide  oxidation  experiments.  In 
these  experiments  the  secondary  flow  consisted  of  N’,  <)•,  and  diethyl  sulfide. 


Condition 

Nominal 

Temperature 

ro 

Flowrate 
of  Main 

N>  flow 
(slpm) 

Flowrate  of 
Secondary 
N-  flow 
(slpm) 

Well- 

mixed 
initial  Oj 
loading 
in 

reaction 

section 

(ppm) 

Well- 

mixed 

initial 

diethyl 

sulfide 

loading 

(ppm) 

Equivalence 

Ratio 

Reynolds 

number 

DOl 

630 

330 

13.3 

733 

100 

1.02 

5020 

DO’ 

670 

341 

14.4 

70S 

100 

1.06 

5060 

D03 

700 

347 

15.2 

695 

100 

1.0S 

5060 

DG4 

740 

350 

16.0 

6X7 

100 

1.09 

4990 

DOS 

700 

347 

16.4 

9042 

100 

0.0X3 

5060 

D06 

740 

350 

13.7 

I960 

100 

0.3S3 

4990 

DO? 

740 

350 

17.3 

9415 

100 

0.0X0 

4990 

Table  6.  Flow  reactor  operating  conditions  for  ethyl  methyl  sulfide  oxidation  experiments. 
In  these  experiments  the  secondary  How  consisted  of  N.,  ():.  and  ethyl  methyl  sulfide. 


Condition 

Nominal 

Temperature 

<°C) 

Flowrate 
of  Mam 
N:  flow 
(slpm) 

Well- 

mixed 

initial 

°2 

loading 

in 

reaction 

section 

(ppm) 

Well- 

mixed 

initial 

ethyl 

methyl 

sulfide 

loading 

(ppm) 

Equivalence 

Ratio 

Reynolds 

number 

EOl 

700 

347 

16.4 

9042 

150 

0.124 

5060 

E02 

740 

350 

17.3 

9415 

150 

0.119 

4990 

IV.B.4.  Res  nils:  species  profiles 

Figures  5  through  1 8  show  experimental  and  calculated  species  profiles  for  the  conditions 
described  in  the  previous  section.  Sampling  locations  have  been  converted  to  residence 
times  using  the  bulk  velocity.  As  mentioned  above,  calculations  were  performed  only  for 


14 


the  diethyl  sulfide  runs,  and  make  use  of  two  different  versions  of  the  mechanism. 
Version  1  is  used  in  the  pyrolysis  simulations,  and  version  2  is  used  in  the  oxidation 
simulations.  Lines  represent  calculations,  while  symbols  represent  experimental  results, 
mainly  with  FT1R.  As  can  be  seen  from  the  figures,  most  experimental  data  points  were 
measured  twice,  with  good  repeatability.  In  some  cases,  both  GC/MS  and  FTIR 
measurement  results  arc  available  and  arc  represented  by  different  symbols.  Agreement 
between  the  two  methods  is  good  in  these  cases.  Included  in  these  figures  arc 
experimental  and  computational  profiles  for  all  species  experimentally  measured  at  levels 
above  I  ppm. 


001  0  02  O.M  OOJ  0  0*  00*  0  0? 


Residence  time  (s) 

Fig.  5.  Major  species  (>  10  ppm)  in  pyrolysis  of  diethyl  sulfide.  Lines  represent  calculations 
with  mechanism  version  I;  symbols  represent  experiments.  Hollow  symbols  represent  FTIR 
data,  w  hile  filled  symbols  represent  GC/MS  data.  Operating  conditions  DPI,  DP2,  DP3, 
and  1)1*4. 
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Fiji.  6.  Minor  species  (<10  ppm)  in  pyrolysis  of  diethyl  sulfide.  Lines  represent  calculations 
with  mechanism  version  1:  symbols  represent  experiments.  Hollow  symbols  represent  FT1R 
data,  while  Tilled  symbols  represent  GC/MS  data.  Operating  conditions  DPI,  DP2,  DP3, 
and  DIM. 
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Fig.  7a  Major  species  (>10  ppm)  for  pyrolysis  of  ethyl  methyl  sulfide  at  630  "C  (operating 
condition  EPI).  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 


(■C/MS  data. 


Residence  time  (s) 


Fig.  7b  Minor  species  (<10  ppm)  for  pyrolysis  of  ethyl  methyl  sulfide  at  630  "C  (operating 
condition  EPI).  FTIR  data  only. 
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Fig.  8a  Major  species  <>10  ppm)  for  pyrolysis  of  ethyl  methyl  sulfide  at  700  "C  (operating 
condition  EP2).  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 
GC/MS  data. 
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Fig.  8I>  Minor  species  (<I0  ppm)  for  pyrolysis  of  ethyl  methyl  sulfide  at  700  "C  (operating 
condition  EP2).  FTIR  data  only. 
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Fig.  9a  Major  species  (>10  ppm)  for  pyrolysis  of  ethyl  methyl  sulfide  at  740  "C  (operating 
condition  EP3).  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 
GC/MS  data. 
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Fig.  9b  Minor  species  (<10  ppm)  for  pyrolysis  of  ethyl  methyl  sulfide  at  740  "C  (operating 
condition  EP3).  Hollow  symbols  represent  FTIR  data,  while  filled  symbols  represent 
GC/MS  data. 
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Fig.  10a  Major  species  (>10  ppm)  for  oxidation  of  diethyl  sulfide  at  630  "C  with  low  O, 
loading  (operating  condition  1)01).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  Hollow  symbols  represent  FT1R  data,  while  filled  symbols 
represent  GC/MS  data. 
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Fig.  10b  Minor  species  (<10  ppm)  for  oxidation  of  diethyl  sulfide  at  630  °C  with  low  <): 
loading  (operating  condition  1)01).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  FT1R  data  only. 
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Fig.  I  la  Major  species  (>10  ppm)  for  oxidation  of  diethyl  sulfide  at  670  "C  with  low  O, 
loading  (operating  condition  1)02).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols 
represent  GC/MS  data. 
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Fig.  lib  Minor  species  (<I0  ppm)  for  oxidation  of  diethyl  sulfide  at  670  “C  with  low  O- 
loading  (operating  condition  1)02).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  FTIR  data  only. 
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Fig.  12a  Major  species  (>10  ppm)  for  oxidation  of  diethyl  sulfide  at  700  "C  with  low  O. 
loading  (operating  condition  1)03).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols 
represent  GC/MS  data. 
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Fig.  12b  Minor  species  (<10  ppm)  for  oxidation  of  diethyl  sulfide  at  700  "C  with  low  O- 
loading  (operating  condition  1)03).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  FTIR  data  only. 
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Fig.  13a  Major  species  (>10  ppm)  for  oxidation  of  diethyl  sulfide  at  740  °C  with  low  O, 
loading  (operating  condition  1)04).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols 
represent  GC/MS  data. 
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Fig.  13b  Minor  species  (<10  ppm)  for  oxidation  of  diethyl  sulfide  at  740  "C  with  low  O- 
loading  (operating  condition  1)04).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  FTIR  data  only. 
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Fig.  14a  Major  species  (>1()  ppm)  for  oxidation  of  diethyl  sulfide  at  700  "C  with  high  O, 
loading  (operating  condition  1)05).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols 
represent  GC/MS  data. 
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Fig.  14b  Minor  species  (<I0  ppm)  for  oxidation  of  diethyl  sulfide  at  700  "C  with  high  O- 
loading  (operating  condition  1)05).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  FTIR  data  only. 
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Fig.  1 5I>  Minor  species  (<10  ppm)  for  oxidation  of  diethyl  sulfide  at  740  "C  with  medium  <)» 
loading  (operating  condition  1)06).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  FTIR  data  only. 


25 


0.01  0.02 

0.03  0.04  0.05  0.06 

Residence  time  (s) 

Fig.  16a  Major  species  (>10  ppm)  for  oxidation  of  diethyl  sulfide  at  740  "C  with  high  O, 
loading  (operating  condition  1)07).  Lines  represent  calculations  with  mechanism  version  2; 
symbols  represent  experiments.  Hollow  symbols  represent  FTIR  data,  while  filled  symbols 
represent  GC/MS  data. 
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Fig.  16b  Minor  species  (<10  ppm  measured)  for  oxidation  of  diethyl  sulfide  at  740  *C  with 
high  <):  loading  (operating  condition  D()7).  Lines  represent  calculations  with  mechanism 
version  2;  symbols  represent  experiments.  FTIR  data  only. 
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Fig.  17b  Minor  species  (<I0  ppm)  for  oxidation  of  ethyl  methyl  sulfide  at  700  *C  with  high 
();  loading  (operating  condition  EOl).  FTIR  data  only 
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Fig.  I8b  Minor  species  (<I0  ppm)  for  oxidation  of  ethyl  methyl  sulfide  at  740  *C  (operating 
condition  EQ2).  FT1R  data  only. 
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/V.B.5.  Results:  other  species 

Sulfur  balances  obtained  from  experimentally  measured  mole  fractions  were  poor;  carbon 
and  hydrogen  balances  were  better,  but  nonetheless  imply  the  presence  significant 
undetected  species.  Figure  19  shows  the  element  balances  for  the  diethyl  sulfide 
pyrolysis  data  set  as  functions  of  destruction  efficiency;  other  data  sets  gave  similar 
results. 


Fig.  19  Element  balances  for  diethyl  sulfide  pyrolysis  data  sets  (operating  conditions  DPI, 
DP2,  DP3.  and  DP4.) 

The  most  likely  cause  of  the  element  balances  is  the  presence  of  major  species  that  arc 
either  undetectable  or  unable  to  pass  through  the  sampling  system.  The  kinetics 
calculations  indicate  the  presence  of  three  major  species  that  arc  not  detected:  hydrogen 
(H:).  thioformaldehydc  (CH’=S),  and  thioacetaldehyde  (CHsCH=S).  Of  these,  hydrogen 
is  not  detectable  by  either  FTIR  or  GC/MS  with  the  current  analytical  method.  The  two 
thioaldehydcs,  on  the  other  hand,  should  be  detectable  by  FTIR  at  the  levels  predicted. 
The  absence  of  their  spectral  features  in  the  FTIR  absorption  spectra  suggests  that  they 
arc  lost  in  the  sampling  system.  One  possibility  was  that  trimers  would  form  from  these 
compounds  [5]  and  would  not  be  volatile  enough  to  travel  through  the  sampling  system; 
this  possibility  was  investigated  as  described  in  the  following  paragraph. 

In  addition,  a  number  of  species  were  observed  experimentally,  by  GC/MS,  but  not 
quantified.  Of  these,  SO’  appears  to  be  present  in  significant  quantities  in  the  oxidation 
runs  and  we  have  plans  to  quantify  it.  Other  species,  for  example:  ethanethiol,  methyl 
thiirane,  ethyl  methyl  disulfide,  and  diethyl  disulfide,  arc  estimated  to  be  present  at  levels 
below  1  ppm  and  do  not  have  a  significant  impact  on  the  element  balances. 

Light  yellow  condensates  were  observed  in  the  sampling  probes  in  the  pyrolysis 
experiments,  and  their  composition  was  investigated  as  follows:  Methylene  chloride 
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(CH’Ch)  was  used  to  rinse  these  sampling  probes,  and  part  of  the  light  yellow  condensate 
was  dissolved  in  the  solvent.  GC/MS  analysis  was  performed  with  a  J&W  DB5  column. 
Although  the  rinsing  and  analysis  method  was  able  to  detect  thioaldehyde  trimer  from  a 
purchased  sample,  no  such  trimer  was  detected  in  the  sample  from  the  probes.  Instead  a 
cyclic  elemental  sulfur  compound,  S*.  was  identified  in  the  rinsatc.  Thus  the  analysis  of 
the  rinsate  provided  no  direct  evidence  of  the  thioaldchydcs  predicted  by  the  kinetic 
calculations.  Clearly,  the  presence  of  sulfur  in  the  wall  deposits  helps  to  explain  the  poor 
sulfur  balance,  but  the  magnitude  of  the  contribution  is  unknown,  as  is  the  route  by  which 
the  cyclic  sulfur  compound  is  formed. 

71  ( 'omggrison  between  calculatigns  and  experiment 

Two  different  versions  of  the  kinetic  mechanism  were  compared  to  predictions.  Version 
1  was  used  for  the  diethyl  sulfide  pyrolysis  calculations,  and  this  version  was  successful 
in  predicting  diethyl  sulfide  destruction.  The  only  major  product  observed 
experimentally  was  ethylene,  which  appears  at  roughly  the  rate  that  the  diethyl  sulfide 
disappears.  The  calculations  underpredict  observed  ethylene  levels  by  approximately 
30%.  Among  the  product  species  observed  and  predicted,  methane  is  overpredicted  by 
approximately  50  to  100%.  Carbon  disulfide,  thiophene,  and  acetylene  arc 
underpredicted.  Ethane  is  underprcdictcd  at  low  temperatures,  but  overpredicted  at  high 
temperatures. 

Version  2  of  the  kinetic  mechanism  was  used  for  the  diethyl  sulfide  oxidation 
calculations.  This  version  underpredicted  destruction  of  diethyl  sulfide  under  all 
conditions.  Ethylene  levels  were  approximately  correct  for  conditions  with  low 
destruction  efficiencies,  and  significantly  underpredicted  under  other  conditions.  CO 
levels  were  overpredicted  by  a  large  factor  under  conditions  of  low  destruction,  when  CO 
levels  were  low.  Agreement  was  better  at  higher  destruction  efficiency.  Minor  species 
C’H6  and  CH4  were  consistently  overpredicted.  Carbon  disulfide  is  underprcdictcd. 

IV.B.7.  Destruction  pathways 

Reaction  pathway  analysis  was  performed  to  identify  the  major  reactions  responsible  for 
the  creation  and  destruction  of  important  species.  The  results  arc  shown  in  Figures  20 
and  2 1  for  diethyl  sulfide  pyrolysis  and  oxidation,  respectively. 


Fig.  20  Major  reaction  pathways  for  diethyl  sulfide  pyrolysis.  Dotted  lines  indicate  radical 
attack  reactions.  Thicker  lines  represent  reactions  with  larger  rates. 
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Fiji.  21  Major  reaction  pathways  for  diethyl  sulfide  oxidation.  Dotted  lines  indicate  radical 
attack  reactions.  Thicker  lines  represent  reactions  with  larger  rates. 


As  the  figures  show,  the  destruction  of  diethyl  sulfide  proceeds  through  basically  the 
same  routes  under  both  oxidative  and  pyrolytic  conditions.  Hydrogen  abstraction  by 
radicals  creates  one  of  two  radicals,  which  then  undergo  beta  scission.  Alternatively, 
unimolecular  dissociation  can  sever  one  of  the  C-S  bonds  of  the  parent  molecule.  The 
products  of  the  dissociation  reactions  undergo  beta  scission. 


V.  REFERENCES 


[1]  M.K.,  Denison.  C.J.  Montgomery.  A.F.  Sarofim.  M.J.  Bockclic,  A.G.  Webster.  R. 
Mellon.  ‘'Advanced  Computational  Modeling  of  Military  Incinerators."  presented  at  the 
Twenty-first  International  Conference  on  Incineration  and  Themial  Treatment 
Technologies.  New  Orleans,  LA.  May  13-17,  2002 

[2]  Sheng.  C.  Ph.D.  Dissertation.  Department  of  Chemical  Engineering.  New  Jersey 
Institute  of  Technology,  2002. 

[3J  X.  Zheng.  E.M.  Fisher.  F.C.  Gouldin.  L.  Zhu.  and  J.W.  Bozzclli,  2008, 

"Experimental  and  Computational  Study  of  Diethyl  Sulfide  Pyrolysis  and  Mechanism." 
accepted  for  publication.  Proceedings  of  the  Combustion  Institute ,  v.  32. 

[4]  Reaction  Design.  2004.  Chcmkin  User  Interface  4.0.  available  at 
hup:  ■.'www.reactiondesi.en.com. 

[5]  Wee  Shong  Chin,  Bee  Wai  Ek,  Chup  Yew  Mok,  Using  Hua  Huang.  J.  Chcm.  Soc. 
Perkin  Trans  2,(1994). 


31 


